Convergent genomic signatures associated with vertebrate viviparity.

BACKGROUND: Viviparity-live birth-is a complex and innovative mode of reproduction that has evolved repeatedly across the vertebrate Tree of Life. Viviparous species exhibit remarkable levels of reproductive diversity, both in the amount of care provided by the parent during gestation, and the ways in which that care is delivered. The genetic basis of viviparity has garnered increasing interest over recent years; however, such studies are often undertaken on small evolutionary timelines, and thus are not able to address changes occurring on a broader scale. Using whole genome data, we investigated the molecular basis of this innovation across the diversity of vertebrates to answer a long held question in evolutionary biology: is the evolution of convergent traits driven by convergent genomic changes? RESULTS: We reveal convergent changes in protein family sizes, protein-coding regions, introns, and untranslated regions (UTRs) in a number of distantly related viviparous lineages. Specifically, we identify 15 protein families showing evidence of contraction or expansion associated with viviparity. We additionally identify elevated substitution rates in both coding and noncoding sequences in several viviparous lineages. However, we did not find any convergent changes-be it at the nucleotide or protein level-common to all viviparous lineages. CONCLUSIONS: Our results highlight the value of macroevolutionary comparative genomics in determining the genomic basis of complex evolutionary transitions. While we identify a number of convergent genomic changes that may be associated with the evolution of viviparity in vertebrates, there does not appear to be a convergent molecular signature shared by all viviparous vertebrates. Ultimately, our findings indicate that a complex trait such as viviparity likely evolves with changes occurring in multiple different pathways.

The genetic basis of a recent transition to live-bearing in marine snails.

Key innovations are fundamental to biological diversification, but their genetic basis is poorly understood. A recent transition from egg-laying to live-bearing in marine snails (Littorina spp.) provides the opportunity to study the genetic architecture of an innovation that has evolved repeatedly across animals. Individuals do not cluster by reproductive mode in a genome-wide phylogeny, but local genealogical analysis revealed numerous small genomic regions where all live-bearers carry the same core haplotype. Candidate regions show evidence for live-bearer-specific positive selection and are enriched for genes that are differentially expressed between egg-laying and live-bearing reproductive systems. Ages of selective sweeps suggest that live-bearer-specific alleles accumulated over more than 200,000 generations. Our results suggest that new functions evolve through the recruitment of many alleles rather than in a single evolutionary step.

Ecological factors and parity mode correlate with genome size variation in squamate reptiles.

BACKGROUND: Evidence of correlation between genome size, the nuclear haploid DNA content of a cell, environmental factors and life-history traits have been reported in many animal species. Genome size, however, spans over three orders of magnitude across taxa and such a correlation does not seem to follow a universal pattern. In squamate reptiles, the second most species-rich order of vertebrates, there are currently no studies investigating drivers of genome size variability. We run a series of phylogenetic generalized least-squares models on 227 species of squamates to test for possible relationships between genome size and ecological factors including latitudinal distribution, bioclimatic variables and microhabitat use. We also tested whether genome size variation can be associated with parity mode, a highly variable life history trait in this order of reptiles. RESULTS: The best-fitting model showed that the interaction between microhabitat use and parity mode mainly accounted for genome size variation. Larger genome sizes were found in live-bearing species that live in rock/sand ecosystems and in egg-laying arboreal taxa. On the other hand, smaller genomes were found in fossorial live-bearing species. CONCLUSIONS: Environmental factors and species parity mode appear to be among the main parameters explaining genome size variation in squamates. Our results suggest that genome size may favour adaptation of some species to certain environments or could otherwise result from the interaction between environmental factors and parity mode. Integration of genome size and genome sequencing data could help understand the role of differential genome content in the evolutionary process of genome size variation in squamates.

Genomic Signatures Associated with Transitions to Viviparity in Cyprinodontiformes.

The transition from oviparity to viviparity has occurred independently over 150 times across vertebrates, presenting one of the most compelling cases of phenotypic convergence. However, whether the repeated, independent evolution of viviparity is driven by redeployment of similar genetic mechanisms and whether these leave a common signature in genomic divergence remains largely unknown. Although recent investigations into the evolution of viviparity have demonstrated striking similarity among the genes and molecular pathways involved across disparate vertebrate groups, quantitative tests for genome-wide convergent have provided ambivalent answers. Here, we investigate the potential role of molecular convergence during independent transitions to viviparity across an order of ray-finned freshwater fish (Cyprinodontiformes). We assembled de novo genomes and utilized publicly available genomes of viviparous and oviparous species to test for molecular convergence across both coding and noncoding regions. We found no evidence for an excess of molecular convergence in amino acid substitutions and in rates of sequence divergence, implying independent genetic changes are associated with these transitions. However, both statistical power and biological confounds could constrain our ability to detect significant correlated evolution. We therefore identified candidate genes with potential signatures of molecular convergence in viviparous Cyprinodontiformes lineages. Motif enrichment and gene ontology analyses suggest transcriptional changes associated with early morphogenesis, brain development, and immunity occurred alongside the evolution of viviparity. Overall, however, our findings indicate that independent transitions to viviparity in these fish are not strongly associated with an excess of molecular convergence, but a few genes show convincing evidence of convergent evolution.

Different Genes are Recruited During Convergent Evolution of Pregnancy and the Placenta.

The repeated evolution of the same traits in distantly related groups (convergent evolution) raises a key question in evolutionary biology: do the same genes underpin convergent phenotypes? Here, we explore one such trait, viviparity (live birth), which, qualitative studies suggest, may indeed have evolved via genetic convergence. There are >150 independent origins of live birth in vertebrates, providing a uniquely powerful system to test the mechanisms underpinning convergence in morphology, physiology, and/or gene recruitment during pregnancy. We compared transcriptomic data from eight vertebrates (lizards, mammals, sharks) that gestate embryos within the uterus. Since many previous studies detected qualitative similarities in gene use during independent origins of pregnancy, we expected to find significant overlap in gene use in viviparous taxa. However, we found no more overlap in uterine gene expression associated with viviparity than we would expect by chance alone. Each viviparous lineage exhibits the same core set of uterine physiological functions. Yet, contrary to prevailing assumptions about this trait, we find that none of the same genes are differentially expressed in all viviparous lineages, or even in all viviparous amniote lineages. Therefore, across distantly related vertebrates, different genes have been recruited to support the morphological and physiological changes required for successful pregnancy. We conclude that redundancies in gene function have enabled the repeated evolution of viviparity through recruitment of different genes from genomic "toolboxes", which are uniquely constrained by the ancestries of each lineage.

Understanding the evolution of viviparity using intraspecific variation in reproductive mode and transitional forms of pregnancy.

How innovations such as vision, flight and pregnancy evolve is a central question in evolutionary biology. Examination of transitional (intermediate) forms of these traits can help address this question, but these intermediate phenotypes are very rare in extant species. Here we explore the biology and evolution of transitional forms of pregnancy that are midway between the ancestral state of oviparity (egg-laying) and the derived state, viviparity (live birth). Transitional forms of pregnancy occur in only three vertebrates, all of which are lizard species that also display intraspecific variation in reproductive phenotype. In these lizards (Lerista bougainvillii, Saiphos equalis, and Zootoca vivipara), geographic variation of three reproductive forms occurs within a single species: oviparity, viviparity, and a transitional form of pregnancy. This phenomenon offers the valuable prospect of watching 'evolution in action'. In these species, it is possible to conduct comparative research using different reproductive forms that are not confounded by speciation, and are of relatively recent origin. We identify major proximate and ultimate questions that can be addressed in these species, and the genetic and genomic tools that can help us understand how transitional forms of pregnancy are produced, despite predicted fitness costs. We argue that these taxa represent an excellent prospect for understanding the major evolutionary shift between egg-laying and live birth, which is a fundamental innovation in the history of animals.

Independent evolutionary transitions to pueriparity across multiple timescales in the viviparous genus Salamandra.

The ability to bear live offspring, viviparity, has evolved multiple times across the tree of life and is a remarkable adaptation with profound life-history and ecological implications. Within amphibians the ancestral reproductive mode is oviparity followed by a larval life stage, but viviparity has evolved independently in all three amphibian orders. Two types of viviparous reproduction can be distinguished in amphibians; larviparity and pueriparity. Larviparous amphibians deliver larvae into nearby ponds and streams, while pueriparous amphibians deliver fully developed juveniles and thus do not require waterbodies for reproduction. Among amphibians, the salamander genus Salamandra is remarkable as it exhibits both inter- and intraspecific variation in the occurrence of larviparity and pueriparity. While the evolutionary relationships among Salamandra lineages have been the focus of several recent studies, our understanding of how often and when transitions between modes occurred is still incomplete. Furthermore, in species with intraspecific variation, the reproductive mode of a given population can only be confirmed by direct observation of births and thus the prevalence of pueriparous populations is also incompletely documented. We used sequence capture to obtain 1,326 loci from 94 individuals from across the geographic range of the genus, focusing on potential reproductive mode transition zones. We also report additional direct observations of pueriparous births for 20 new locations and multiple lineages. We identify at least five independent transitions from the ancestral mode of larviparity to pueriparity among and within species, occurring at different evolutionary timescales ranging from the Pliocene to the Holocene. Four of these transitions occurred within species. Based on a distinct set of markers and analyses, we also confirm previous findings of introgression between species and the need for taxonomic revisions in the genus. We discuss the implications of our findings with respect to the evolution of this complex trait, and the potential of using five independent convergent transitions for further studies on the ecological context in which pueriparity evolves and the genetic architecture of this specialized reproductive mode.

The Evolution of Viviparity in Vertebrates.

In the vertebrate tree of life, viviparity or live birth has independently evolved many times, resulting in a rich diversity of reproductive strategies. Viviparity is believed to be a mode of reproduction that evolved from the ancestral condition of oviparity or egg laying, where most of the fetal development occurs outside the body. Today, there is not a simple model of parity transition to explain this species-specific divergence in modes of reproduction. Most evidence points to a gradual series of evolutionary adaptations that account for this phenomenon of reproduction, elegantly displayed by various viviparous squamates that exhibit placentae formed by the appositions of maternal and embryonic tissues, which share significant homology with the tissues that form the placenta in therian mammals. In an era where the genomes of many vertebrate species are becoming available, studies are now exploring the molecular basis of this transition from oviparity to viviparity, and in some rare instances its possible reversibility, such as the Australian three-toed skink (Saiphos equalis). In contrast to the parity diversity in squamates, mammals are viviparous with the notable exception of the egg-laying monotremes. Advancing computational tools coupled with increasing genome availability across species that utilize different reproductive strategies promise to reveal the molecular underpinnings of the ancestral transition of oviparity to viviparity. As a result, the dramatic changes in reproductive physiology and anatomy that accompany these parity changes can be reinterpreted. This chapter will briefly explore the vertebrate modes of reproduction using a phylogenetic framework and where possible highlight the role of potential candidate genes that may help explain the polygenic origins of live birth.

Parallel Genomic Changes Drive Repeated Evolution of Placentas in Live-Bearing Fish.

The evolutionary origin of complex organs challenges empirical study because most organs evolved hundreds of millions of years ago. The placenta of live-bearing fish in the family Poeciliidae represents a unique opportunity to study the evolutionary origin of complex organs, because in this family a placenta evolved at least nine times independently. It is currently unknown whether this repeated evolution is accompanied by similar, repeated, genomic changes in placental species. Here, we compare whole genomes of 26 poeciliid species representing six out of nine independent origins of placentation. Evolutionary rate analysis revealed that the evolution of the placenta coincides with convergent shifts in the evolutionary rate of 78 protein-coding genes, mainly observed in transporter- and vesicle-located genes. Furthermore, differences in sequence conservation showed that placental evolution coincided with similar changes in 76 noncoding regulatory elements, occurring primarily around genes that regulate development. The unexpected high occurrence of GATA simple repeats in the regulatory elements suggests an important function for GATA repeats in developmental gene regulation. The distinction in molecular evolution observed, with protein-coding parallel changes more often found in metabolic and structural pathways, compared with regulatory change more frequently found in developmental pathways, offers a compelling model for complex trait evolution in general: changing the regulation of otherwise highly conserved developmental genes may allow for the evolution of complex traits.

Evolution of reproductive strategies in the species-rich land snail subfamily Phaedusinae (Stylommatophora: Clausiliidae).

Most of the present knowledge on animal reproductive mode evolution, and possible factors driving transitions between oviparity and viviparity is based on studies on vertebrates. The species rich door snail (Clausiliidae) subfamily Phaedusinae represents a suitable and unique model for further examining parity evolution, as three different strategies, oviparity, viviparity, and the intermediate mode of embryo-retention, occur in this group. The present study reconstructs the evolution of reproductive strategies in Phaedusinae based on time-calibrated molecular phylogenetics, reproductive mode examinations and ancestral state reconstruction. Our phylogenetic analysis employing multiple mitochondrial and nuclear markers identified a well-supported clade (including the tribes Phaedusini and Serrulinini) that contains species exhibiting various reproductive strategies. This clade evolved from an oviparous most recent common ancestor according to our reconstruction. All non-oviparous taxa are confined to a highly supported subclade, coinciding with the tribe Phaedusini. Both oviparity and viviparity occur frequently in different lineages of this subclade that are not closely related. During Phaedusini diversification, multiple transitions in reproductive strategy must have taken place, which could have been promoted by a high fitness of embryo-retaining species. The evolutionary success of this group might result from the maintenance of various strategies.

Chromosome-Level Assembly of the Common Lizard (Zootoca vivipara) Genome.

Squamate reptiles exhibit high variation in their phenotypic traits and geographical distributions and are therefore fascinating taxa for evolutionary and ecological research. However, genomic resources are very limited for this group of species, consequently inhibiting research efforts. To address this gap, we assembled a high-quality genome of the common lizard, Zootoca vivipara (Lacertidae), using a combination of high coverage Illumina (shotgun and mate-pair) and PacBio sequencing data, coupled with RNAseq data and genetic linkage map generation. The 1.46-Gb genome assembly has a scaffold N50 of 11.52 Mb with N50 contig size of 220.4 kb and only 2.96% gaps. A BUSCO analysis indicates that 97.7% of the single-copy Tetrapoda orthologs were recovered in the assembly. In total, 19,829 gene models were annotated to the genome using a combination of ab initio and homology-based methods. To improve the chromosome-level assembly, we generated a high-density linkage map from wild-caught families and developed a novel analytical pipeline to accommodate multiple paternity and unknown father genotypes. We successfully anchored and oriented almost 90% of the genome on 19 linkage groups. This annotated and oriented chromosome-level reference genome represents a valuable resource to facilitate evolutionary studies in squamate reptiles.

Viviparous Reptile Regarded to Have Temperature-Dependent Sex Determination Has Old XY Chromosomes.

The water skinks Eulamprus tympanum and Eulamprus heatwolei show thermally induced sex determination where elevated temperatures give rise to male offspring. Paradoxically, Eulamprus species reproduce in temperatures of 12-15 °C making them outliers when compared with reptiles that use temperature as a cue for sex determination. Moreover, these two species are among the very few viviparous reptiles reported to have thermally induced sex determination. Thus, we tested whether these skinks possess undetected sex chromosomes with thermal override. We produced transcriptome and genome data for E. heatwolei. We found that E. heatwolei presents XY chromosomes that include 14 gametologs with regulatory functions. The Y chromosomal region is 79-116 Myr old and shared between water and spotted skinks. Our work provides clear evidence that climate could be useful to predict the type of sex determination systems in reptiles and it also indicates that viviparity is strictly associated with sex chromosomes.

Molecular mechanisms underlying milk production and viviparity in the cockroach, Diploptera punctata.

Viviparous reproduction is characterized by maternal retention of developing offspring within the reproductive tract during gestation, culminating in live birth. In some cases, a mother will provide nutrition beyond that present in the yolk; this is known as matrotrophic viviparity. While this phenomenon is best associated with mammals, it is observed in insects such as the viviparous cockroach, Diploptera punctata. Female D. punctata carry developing embryos in the brood sac, a reproductive organ that acts as both a uterus and a placenta by protecting and providing a nutritive secretion to the intrauterine developing progeny. While the basic physiology of D. punctata pregnancy has been characterized, little is known about the molecular mechanisms underlying this phenomenon. This study combined RNA-seq analysis, RNA interference, and other assays to characterize molecular and physiological changes associated with D. punctata reproduction. A comparison of four stages of the female reproductive cycle and males revealed unique gene expression profiles corresponding to each stage and between sexes. Differentially regulated transcripts of interest include the previously identified family of milk proteins and transcripts associated with juvenile hormone metabolism. RNA interference and methoprene application experiments established the potential impacts of bothbreakdown and synthesis reduction of juvenile hormone in maintaining pregnancy in D. punctata. These studies provide the comprehensive molecular mechanisms associated with cockroach viviparity, which will be a critical resource for comparative purposes among viviparity in insect systems.

Emergence of an evolutionary innovation: Gene expression differences associated with the transition between oviparity and viviparity.

Our understanding of the evolution of complex biological traits is greatly advanced by examining taxa with intermediate phenotypes. The transition from oviparity (egg-laying) to viviparity (live-bearing) has occurred independently in many animal lineages, but there are few phenotypic intermediates. The lizard Saiphos equalis exhibits bimodal reproduction, with some viviparous populations, and other oviparous populations with long egg-retention, a rare trait where most of embryonic development occurs inside the mother prior to late ovipositioning. We posit that oviparous S. equalis represent an intermediate form between "true" oviparity and viviparity. We used transcriptomics to compare uterine gene expression in these two phenotypes, and provide a molecular model for the genetic control and evolution of reproductive mode. Many genes are differentially expressed throughout the reproductive cycle of both phenotypes, which have clearly different gene expression profiles overall. The differentially expressed genes within oviparous and viviparous individuals have broadly similar biological functions putatively important for sustaining embryos, including uterine remodelling, respiratory gas and water exchange, and immune regulation. These functional similarities indicate either that long egg-retention is an exaptation for viviparity, or might reflect parallel evolution of similar gravidity-related changes in gene expression in long egg-retention oviparity. In contrast, gene expression changes across the reproductive cycle of long egg-retaining oviparous S. equalis are dramatically different from those of "true" oviparous skinks (such as Lampropholis guichenoti), supporting our assertion that oviparous S. equalis exhibit an intermediate phenotype between "true" oviparity and viviparity.

The Genomes of the Livebearing Fish Species Poeciliopsis retropinna and Poeciliopsis turrubarensis Reflect Their Different Reproductive Strategies.

The evolution of a placenta is predicted to be accompanied by rapid evolution of genes involved in processes that regulate mother-offspring interactions during pregnancy, such as placenta formation, embryonic development, and nutrient transfer to offspring. However, these predictions have only been tested in mammalian species, where only a single instance of placenta evolution has occurred. In this light, the genus Poeciliopsis is a particularly interesting model for placenta evolution, because in this genus a placenta has evolved independently from the mammalian placenta. Here, we present and compare genome assemblies of two species of the livebearing fish genus Poeciliopsis (family Poeciliidae) that differ in their reproductive strategy: Poeciliopsis retropinna which has a well-developed complex placenta and P. turrubarensis which lacks a placenta. We applied different assembly strategies for each species: PacBio sequencing for P. retropinna (622-Mb assembly, scaffold N50 of 21.6 Mb) and 10× Genomics Chromium technology for P. turrubarensis (597-Mb assembly, scaffold N50 of 4.2 Mb). Using the high contiguity of these genome assemblies and near-completeness of gene annotations to our advantage, we searched for gene duplications and performed a genome-wide scan for genes evolving under positive selection. We find rapid evolution in major parts of several molecular pathways involved in parent-offspring interaction in P. retropinna, both in the form of gene duplications as well as positive selection. We conclude that the evolution of the placenta in the genus Poeciliopsis is accompanied by rapid evolution of genes involved in similar genomic pathways as found in mammals.

Sequencing of the black rockfish chromosomal genome provides insight into sperm storage in the female ovary.

Black rockfish (Sebastes schlegelii) is an economically important viviparous marine teleost in Japan, Korea, and China. It is characterized by internal fertilization, long-term sperm storage in the female ovary, and a high abortion rate. For better understanding the mechanism of fertilization and gestation, it is essential to establish a reference genome for viviparous teleosts. Herein, we used a combination of Pacific Biosciences sequel, Illumina sequencing platforms, 10× Genomics, and Hi-C technology to obtain a genome assembly size of 848.31 Mb comprising 24 chromosomes, and contig and scaffold N50 lengths of 2.96 and 35.63 Mb, respectively. We predicted 39.98% repetitive elements, and 26,979 protein-coding genes. S. schlegelii diverged from Gasterosteus aculeatus ∼32.1-56.8 million years ago. Furthermore, sperm remained viable within the ovary for up to 6 months. The glucose transporter SLC2 showed significantly positive genomic selection, and carbohydrate metabolism-related KEGG pathways were significantly up-regulated in ovaries after copulation. In vitro suppression of glycolysis with sodium iodoacetate reduced sperm longevity significantly. The results indicated the importance of carbohydrates in maintaining sperm survivability. Decoding the S. schlegelii genome not only provides new insights into sperm storage; additionally, it is highly valuable for marine researchers and reproduction biologists.

Comparative landscape genetics reveals the evolution of viviparity reduces genetic connectivity in fire salamanders.

Evolutionary changes in reproductive mode may affect co-evolving traits, such as dispersal, although this subject remains largely underexplored. The shift from aquatic oviparous or larviparous reproduction to terrestrial viviparous reproduction in some amphibians entails skipping the aquatic larval stage and, thus, greater independence from water. Accordingly, amphibians exhibiting terrestrial viviparous reproduction may potentially disperse across a wider variety of suboptimal habitats and increase population connectivity in fragmented landscapes compared to aquatic-breeding species. We investigated this hypothesis in the fire salamander (Salamandra salamandra), which exhibits both aquatic- (larviparity) and terrestrial-breeding (viviparity) strategies. We genotyped 426 larviparous and 360 viviparous adult salamanders for 13 microsatellite loci and sequenced a mitochondrial marker for 133 larviparous and 119 viviparous individuals to compare population connectivity and landscape resistance to gene flow within a landscape genetics framework. Contrary to our predictions, viviparous populations exhibited greater differentiation and reduced genetic connectivity compared to larviparous populations. Landscape genetic analyses indicate viviparity may be partially responsible for this pattern, as water courses comprised a significant barrier only in viviparous salamanders, probably due to their fully terrestrial life cycle. Agricultural areas and, to a lesser extent, topography also decreased genetic connectivity in both larviparous and viviparous populations. This study is one of very few to explicitly demonstrate the evolution of a derived reproductive mode affects patterns of genetic connectivity. Our findings open avenues for future research to better understand the eco-evolutionary implications underlying the emergence of terrestrial reproduction in amphibians.

Genomic and transcriptomic investigations of the evolutionary transition from oviparity to viviparity.

Viviparous (live-bearing) vertebrates have evolved repeatedly within otherwise oviparous (egg-laying) clades. Over two-thirds of these changes in vertebrate reproductive parity mode happened in squamate reptiles, where the transition has happened between 98 and 129 times. The transition from oviparity to viviparity requires numerous physiological, morphological, and immunological changes to the female reproductive tract, including eggshell reduction, delayed oviposition, placental development for supply of water and nutrition to the embryo by the mother, enhanced gas exchange, and suppression of maternal immune rejection of the embryo. We performed genomic and transcriptomic analyses of a closely related oviparous-viviparous pair of lizards (Phrynocephalus przewalskii and Phrynocephalus vlangalii) to examine these transitions. Expression patterns of maternal oviduct through reproductive development of the egg and embryo differ markedly between the two species. We found changes in expression patterns of appropriate genes that account for each of the major aspects of the oviparity to viviparity transition. In addition, we compared the gene sequences in transcriptomes of four oviparous-viviparous pairs of lizards in different genera (Phrynocephalus, Eremias, Scincella, and Sphenomorphus) to look for possible gene convergence at the sequence level. We discovered low levels of convergence in both amino acid replacement and evolutionary rate shift. This suggests that most of the changes that produce the oviparity-viviparity transition are changes in gene expression, so occasional reversals to oviparity from viviparity may not be as difficult to achieve as has been previously suggested.

Conservation of Sex-Linked Markers among Conspecific Populations of a Viviparous Skink, Niveoscincus ocellatus, Exhibiting Genetic and Temperature-Dependent Sex Determination.

Sex determination systems are exceptionally diverse and have undergone multiple and independent evolutionary transitions among species, particularly reptiles. However, the mechanisms underlying these transitions have not been established. Here, we tested for differences in sex-linked markers in the only known reptile that is polymorphic for sex determination system, the spotted snow skink, Niveoscincus ocellatus, to quantify the genomic differences that have accompanied this transition. In a highland population, sex is determined genetically, whereas in a lowland population, offspring sex ratio is influenced by temperature. We found a similar number of sex-linked loci in each population, including shared loci, with genotypes consistent with male heterogamety (XY). However, population-specific linkage disequilibrium suggests greater differentiation of sex chromosomes in the highland population. Our results suggest that transitions between sex determination systems can be facilitated by subtle genetic differences.

Uterine and eggshell modifications associated with the evolution of viviparity in South American water snakes (Helicops spp.).

The evolution of viviparity requires eggshell thinning to bring together the maternal uterus and extraembryonic membranes to form placentae for physiological exchanges. Eggshell thinning likely involves reduced activity of the uterine glands that secrete it. We tested these hypotheses by comparing the uterine and eggshell structure and histochemistry among oviparous and viviparous water snakes (Helicops) using phylogenetic methods. Eggshell thinning occurred convergently in all three origins of viviparity in Helicops and was accomplished by the loss of the mineral layer and thinning of the shell membrane. Uterine glands secrete the shell membrane in both oviparous and viviparous Helicops. These glands increase during vitellogenesis regardless of the reproductive mode, but they always reach smaller sizes in viviparous forms. As there is no phylogenetic signal in eggshell thickness and gland dimensions, we conclude that interspecific differences are related to reproductive mode and not phylogeny. Therefore, our results support the hypothesis that eggshell thinning is associated with the evolution of viviparity and that such thinning result from a reduction in gland size in viviparous taxa. Interestingly, the shell membrane thickness of viviparous females of the reproductively bimodal Helicops angulatus is intermediate between their oviparous and viviparous congeners. Thus, although eggshell thinning is required by the evolution of viviparity, a nearly complete loss of this structure is not. However, uterine gland dimensions are similar across viviparous Helicops. Fewer glands or their functional repurposing may explain the thinner shell membrane in viviparous species of Helicops in comparison to viviparous females of the bimodal H. angulatus.